CO adsorption on Cu(211) surface: first principle and STM study. 
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Chemisorption of CO on the stepped Cu(211) surface is studied within ab-initio density func- 
tional theory (DFT) and scanning tunneling microscopy (STM) imaging as well as manipulation 
experiments. Theoretically we focus on the experimentally observed ordered (2x1) and (3x1) CO- 
phases at coverages = |, i and | monolayer (ML). To obtain also information for isolated CO 
molecules found randomly distributed at low coverages, we also performed calculations for a hy- 
pothetical (3x1) phase with O = | ML. The adsorption geometry, the stretching frequencies, the 
work functions and adsorption energies of the CO molecules in the different phases are presented and 
compared to experimental data. Initially and up to a coverage of § ML CO adsorbs upright on the 
on-top sites at step edge atoms. Determining the most favorable adsorption geometry for the | ML 
ordered phase turned out to be nontrivial, both from the experimental and the theoretical point of 
view. Experimentally, both top-bridge and top-top configurations were reported, whereby only the 
top-top arrangement was firmly established. The calculated adsorption energies and the stretching 
frequencies favor the top-bridge configuration. The possible existence of both configurations at | 
ML is critically discussed on the basis of the presently accessible experimental and theoretical data. 
In addition, we present observations of STM manipulation experiments and corresponding theoret- 
ical results, which show that CO adsorbed on-top of a single Cu-adatom, which is manipulated to a 
location close to the lower step edge, is more strongly bound than CO on-top of a step edge atom. 

PACS numbers: 68.43.-h, 71. 15. Mb, 82.37.Gk, 82.65.+r 
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INTRODUCTION 



Experimental and theoretical studies of chemisorption 
on vicinal surfaces are of great current interest due to 
close connection to real catalytic substrates which ex- 
hibit beside other defects many steps. Thanks to scan- 
ning tunneling microscopy (STM) a breakthrough in de- 
termining the morphology of substrates and ovcrlaycrs 
was achieved in recent years Carbon monoxide on 
Cu(211) was studied with respect to ordered overlayers 
and also in connection with the possibility to laterally 
manipulate adsorbates on the surface by STM 0, 0|- 
Meyer et al. 4] found that the Cu(211) surface does 
not reconstruct and that there are several ordered CO 
structures upon adsorption. We concentrate here on the 
(2x1) and the (3x1) phases with respective coverages of 
6=|, i and | ML. STM imaging of disordered adsorp- 
tion at very low coverage and of the (2x1) phase revealed 
that the CO molecules adsorb on top of Cu atoms at 
the step edges. The high coverage (3x1) phase is more 
complicated: upon preparing an adlayer at 100 K, the ad- 
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molecule configuration appears to be top-bridge, whereas 
in small islands prepared at 15 K by lateral manipula- 
tion of the molecules they occupy top-top arrangements 
within the (3x1) unit cell [5j. 

Braun et al. |(j have observed a significant decrease in 
the energy of the frustrated translation vibrations par- 
allel to the surface (T-mode) of isolated CO molecules 
adsorbed on the intrinsic step sites of Cu(211) as com- 
pared to CO on flat Cu(lll) and Cu(OOl) surfaces. They 
propose a tilting of the CO molecules due to the presence 
of these intrinsic steps as a possible explanation for this 
energy decrease. 

There have been several theoretical investigations of 
the geometric structure of clean Cu(211) 0, U and also 
Cu(nll) surfaces [U El Very recently, reviews of elec- 
tronic, structural, vibrational and thermodynamic char- 
acterization of vicinal surfaces were published [UIH E3- 
On the other hand, calculations on the adsorption of CO 
molecules on Cu(211) are less numerous [l4|- Further- 
more, theoretical investigations of vibrational and ener- 
getic properties of the CO adsorption system are scarce 
or missing. Adsorption of CO on Cu(211) is also inter- 
esting from the fundamental point of view as it offers an 
occasion to study adsorbate re orga nization as proposed 
for NO on the Pd(211) surface [13. 

In this work we carried out a systematic theoretical 
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investigation of the clean copper (211) surface and of the 
adsorption of CO molecules. In section 11111 we examine 
the clean Cu(211) surface (multilayer relaxations, work- 
functions). The remainder of the paper is organized as 
follows. Sect ion Hvl is divided into several parts where we 
describe structural, vibrational and energetic properties 
of the CO molecules adsorbed on the Cu(211) surface. 
Moreover, subsection IIV Dl provides calculated STM im- 
ages of the probable structural phases at 6 = § ML. We 
present and discuss our results parallclly to create a pic- 
ture of the trends for clean Cu(211) and CO molecules 
adsorbed on this surface. Finally, section provides a 
short summary of the main results and conclusions. 



II. METHODOLOGY 

In this work, a slab approach is adapted to describe 
the metallic surface. Twelve (resp. 14 layers for clean 
surface) (211) crystallographic layers are repeated pe- 
riodically in a supercell geometry. The vicinal (211) 
surface consists of (111) oriented terrace planes and 
monoatomic (100) steps according to the Somorjai nota- 
tion (3(111) x (100)) 16]. The macroscopic surface plane 
is rotated by 19.5° from the (111) plane. The distance 
between the steps is 6.25 A according to the theoretical 
bulk lattice parameter &c u — 3.664 A of fee Cu. The 
structure of the clean Cu(211) surface is shown in Fig. ^ 
We consider CO molecules in the (2x1) and the (3x1) 
unit cells. 

The calculations presented in this study are performed 
within a plane-wave density functional framework. We 
have used the Vienna Ab-initio Simulation Package 
(VASP) mmi and employed projected augmented-wave 
(PAW) potentials |19l l20l] . For exchange and correlation 
the functional proposed by Perdew and Zunger [21j is 
used, adding (semi-local) generalized gradient approxi- 
mation (GGA) of the PW91 Further, the RPBE 
exchange-correlation functional [23j is used, since it is 
more accurate for the CO adsorption on the Cu(lll) sur- 
face [24|. It is well known that in some cases DFT cal- 
culations fail to produce the correct adsorption site for 
CO on transition or noble metals surfaces [24ll25| . It has 
been shown that this failure arises from the underestima- 
tion of the HOMO-LUMO gap of the CO molecule favor- 
ing backdonation to the antibonding 2n* state and ad- 
sorption in a site with higher coordination j2(| . Very re- 
cently, Kresse et al. [2(j have shown that the correct site- 
preference can be achieved by adding a Hubbard-like on- 
site Coulomb repulsion [13,123] to the molecular Hamilto- 
nian of CO leading to an increased HOMO-LUMO gap. 
We have applied the special GGA+U method [2(| to ex- 
amine the effect of stronger electronic correlation on the 
prediction of the correct adsorption site. 

We have used (4x6x1) and (4x4x1) k-point meshes 
according to Monkhorst-Pack for the (2x1) and (3x1) 
unit cells, respectively. The adsorption energy (E a d s ) 
was defined as the difference between the total energy of 



the studied system (Cu(211) surface with CO) and the 
individual components of the system (Cu(211) surface 
and CO) [!|. 

The atoms in the first five layers and C and O atoms of 
the CO molecule were allowed to relax. The geometry op- 
timizations were caried according to conjugate-gradient 
algorithm. 

The finite-difference method is used for the calculation 
of stretching frequencies. In the finite-difference tech- 
nique, total energies and Hellmann-Fcynman forces are 
evaluated as a function of the atomic displacements from 
their equilibrium positions; C and O atoms are moved 
in the direction of carthesian axis (x,y,z) by 0.04 A and 
the substrate is frozen. The Hessian dynamical matrix 
is created from the second derivatives of the total energy 
with respect to the positions of CO. The matrix is di- 
agonalized and normalized eigenvectors denote modes of 
the vibrations and corresponding eigenvalues are propor- 
tional to stretching frequencies of the vibrational modes 

We have applied a Tersoff-Hamman approach to cal- 
culate STM images from the ab-initio calculations where 
the STM contour plot is approximated by the charge den- 
sity around the Fermi level (± ~ 50 meV) [3lJ . Constant 
current topographs are approximated by constant charge 
density iso-surfaces. The iso-surfaces are given for the 
charge density value of e -5 eV/A 3 . 



III. CLEAN Cu(211) SURFACE 

We have calculated multilayer relaxations of the clean 
Cu(211) surface and also of clean Cu(211) with one Cu- 
adatom on the (111) terrace adjacent to a step (see cross 
in Fig. ^as well as Fig. [31 model a). For the calcula- 
tions, the lateral distance between the extra Cu-adatoms 
was chosen to be 3b (see Fig. ^1 for the extra Cu atoms 
in the direction of the steps. The changes of the in- 
terlayer spacing from both experiment (LEED |32j) and 
theoretical studies are tabulated in Table Q] Notice that 
these changes are not strictly oscillatory. The earlier and 
current theoretical and experimental studies predict the 
same trends for the interlayer relaxation (-,-,+) for four 
surface layers. Such behaviour has been obtained in a 
FLAPW study @, but the authors report a more pro- 
nounced contraction of the first interlayer distance as 
compared to the experimental results and the present 
theoretical study. On the other hand, the LDA Q and 
current GGA calculations give very similar results which 
are moreover close to the experiment. Inclusion of a Cu- 
adatom next to a step on the Cu(211) surface expands 
only its local environment which gives rise to a small 
buckling of the layer below the Cu-adatom. 

The workfunction of clean Cu(211) ($cu(2ii)) was cal- 
culated to be 4.46 eV. It is lower than the workfunction 
for Cu(lll) ($ Ctl ( ln) =5.00 eV) 0], in accordance with 
Smoluchowski's consideration of the smoothening effect 
of the surface charge density that is responsible for the 
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d=6.35A Step plane( 100) 




plane (111) Y (0Tl) 



FIG. 1: Top and side views of vicinal Cu(211) = (3(lll) x (100)). Typical parallelograms for (111) terraces and rectangles for 
(100) steps are indicated. In addition we show (2x1) and (3x1) cells for the CO adsorption structures. 
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LDA+norm" [331 
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-2.1 


-1.6 


+1.5 


-0.1 


LDA [7] 


-14.4 


-10.7 


+10.9 


-3.8 


-2.3 


+1.7 


-1.0 


FPLAPW [8J 


-28.4 


-3.0 


+15.3 


-6.6 


+0.7 


+3.0 


+0.0 


EAM 6 [34] 


-10.3 


-5.4 


+7.3 


-5.7 


-1.2 


+4.0 


-2.6 


EAM [35] 


-10.3 


-5.1 


+7.3 


-5.6 


-1.1 






LEED [32] 


-14.9±4.1 -10.8±4.1 +8.1±4.1 



a Local Density Approximation with with norm-conserving, nonlocal pseudopotentials 
^Embedded-atom method 



TABLE I: Relaxations of the 14 layer clean Cu(211) surface and the clean Cu(211) surface with one Cu-adatom on the terrace 
close to a step (see Fig. |3Ji). The bulk interlayer spacing is dbuik = 0.748 A for Cu(211). Fig. ^indicates the indexing of the 
individual layers. 



workfunction reduction on more open surfaces |36|. Ad- 
ditional Cu atoms on (111) terraces close to the steps 
increase the overall dipole moment and consequently de- 
crease the "workfunction. Summarizing, the interlayer re- 
laxations are well described by our theoretical tool and 
the trend (-,-,+) is consistent with experimental obser- 
vation. 



IV. CO ADSORPTION 

A. Structural phases 

We have searched through the literature to determine 
the most probable adsorption sites for the different cov- 
erages. The majority of investigations claim that CO 
stands more or less upright on the atomic rows forming 
the steps. As credible sites we therefore considered top 
and bridge sites on the upper parts of the step edges. 
We have theoretically examined \ and | ML coverage 



phases of CO molecules as well as a hypothetical = | 
ML phase. We found in the literature also chance for 
CO-adsorption on the (111) surface facet with the CO 
axis almost parallel to the substrate |37||. but we were 
not able to observe a minimum of the potential energy 
surface for this configuration. 



B. Structural and energetic properties 

Our calculations indicate that CO molecules at the 
steps adsorb upright with respect to the macroscopic sur- 
face plane as long as they do not occupy nearest neigh- 
bor sites at the step edges as in one of the possible con- 
figurations in the high coverage (3x1) phase. The CO 
molecules affect only substrate atoms around their local 
environments and reduce the inward relaxation of atoms. 
The same trend has been observed for CO adsorption 
on transition and noble metal surfaces (24[. The sketch 
of the calculated structural and vibrational properties is 
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shown in Fig. [2] The values of the structural properties 
for all calculated adsorption sites are presented in Table 

EH 

The CO bond length increases with the coordination 
due to larger occupation of the 2ir* orbital [24j. There 
are only small changes of the CO bond length with in- 
creasing coverage. Furthermore, the CO molecule ap- 
proaches closer the surface as the coordination increases. 
The corresponding Cu-C bond length increases with the 
coordination. 

The bonds that the CO molecules form with Cu atoms 
reduce the Cu-Cu binding, which is the major reason 
for the buckling of the first substrate layer. The buck- 
ling quickly disappears and is almost zero in the third 
layer. The buckling in the layers influences interlayer re- 
laxations in an averaged way; consequently the expansion 
of the first layer is similar for the bridge site at -| ML 
and the top site at g ML coverages. The reduction of 
the second layer is similar for all coverages (except CO 
adsorption on Cu-adatom). 

The workfunction of clean Cu(211) is 4.46 eV and in- 
creases upon filling top and bridge sites with CO. More- 
over, the workfunction increases with CO coverage and 
also with CO coordination. For the bridge sites the in- 
crease in workfunction is by ~ 0.2 eV larger than for top 
sites in the (3x1) and by ~ 0.4 eV for the (2x1) phase. 
Additionaly, the workfunction for top-top configuration 
of NO in (3x1) cell is 0.18 eV lower than the workfunc- 
tion for top-bridge configuration. This increase may be 
of further help in determining the optimal adsorption site 
by respective experiments. 

To determine the most favorable structures from a cho- 
sen set of sites at a given coverage we present the adsorp- 
tion energies in Table ITTT1 The difference in the adsorp- 
tion energy between top and bridge site at O = | ML 
and i ML is ~ 20 meV for PW91 exchange-correlation 
functional, which is small to distinguish the optimal site. 
We observe, for the first time, that the RPBE exchange- 
correlation functional changes the site preference and 
gives a different answer concerning the site preference as 
PW91. The difference in the RPBE adsorption energy 
between top and bridge site is not large (~ 30 meV), 
but assigns the top site as the optimal site for i and \ 
ML coverage. Our conclusions are similar to the earlier 
semi- empirical calculations of Marinica et al. [H]- Ear- 
lier work on CO on the close-packed transition and noble 
metal surfaces led to the conclusion that all exchange- 
correlation functionals give the same site prediction and 
the difference among them is mainly in the absolute value 
of the adsorption energy [24| . 

Underestimation of the one-electron gap between 
HOMO (5fi) and LUMO (27r*) is the major problem in 
the theoretical description of CO molecules. The gap 
can be enlarged by inclusion of an artificial Coulomb re- 
pulsion U in the calculation in the way described in the 
paper of Kresse et al. [2(| [3(| . The adsorption energy 
decreases linearly with the U parameter and the slope 
increases with coordination. We have taken U = 1.0 eV 



that leads to the top site as the favorite adsorption site 
on the Cu(lll) surface ^cJj m accordance with exper- 
iment. In our calculations this again changes the pref- 
erence from the bridge site to the top site for the CO 
adsorption in the hypothetical (3x1) and the (2x1) unit 
cell. The CO molecule in the on-top site binds stronger 
by ~ 130 meV to the substrate. This also indicates 
that PW91 exchange-correlation functional can not only 
greatly overestimate the adsorption energy, but can also 
lead to a qualitatively different result. 

It should be emphasized that the CO molecule on top 
of an Cu-adatom located at the lower step edge binds 
stronger by 0.16 eV for PW91, 0.16 eV for RPBE and 
0.18 eV for GGA+U than CO at top sites on the steps. 
This is in accordance with experimental STM manipula- 
tion results which are outlined in Fig. [31 In Fig. [3J<Y a 
configuration is shown which was created artificially by 
manipulating a Cu-adatom to an intrinsic step edge of 
the Cu(211) as depicted in the model|3Ji. A CO molecule 
located initially on top of an atom at the step was then 
pushed along the step edge. It was found that the CO 
goes on top of the Cu-adatom as depicted in Fig. [3)3 and 
model b. Upon manipulating the CO molecule further 
to move it back to the step edge (see Fig. |2P and model 
c), it was found that a larger force was required. As di- 
rect force measurements with the STM are not possible, 
the force is measured semi-quantitatively by the manipu- 
lation resistivity: Smaller resistivities imply larger forces 
because the tip is closer to the manipulated particle. The 
respective tunneling resistivities were 400 kfi and 270 
k£7. Thus the experimental result indicates an enhanced 
bonding of the CO to the single Cu-adatom as compared 
to a CO molecule sitting on top of a Cu-atom in the row 
of atoms building up the step. This experimental obser- 
vation is strongly corroborated by the present theoretical 
results. This result is also of appreciable practical im- 
portance as it proves, that the CO molecules transfered 
deliberately to the tip apex by vertical manipulation |4l| 
prefer the adsorption at the most exposed site, i.e. the 
single metal atom at the tip apex. 



C. Vibrational properties 

We present the calculated vibrational frequencies of 
CO molecules adsorbed on Cu(211) in Table ifVl The C- 
O stretching frequencies (yc-o) agree with the already 
proposed trend that the stretching frequency strongly de- 
pends on coordination. In our case the stretching fre- 
quency for the on-top site is 2049 cm -1 and 1925 cm -1 
for the bridge site at a coverage of | ML, i.e. red-shifted 
by 87 cm -1 and 211 cm -1 compared to the calculated 
gas-phase value of 2136 cm -1 |24|. 

vc-o increases with the coverage for the on-top sites 
and decreases with increased coordination (bridge sites). 
A CO molecule adsorbed on-top of the Cu single adatom 
adjacent to a step edge vibrates in a steep potential as 
demonstrated by the highest C-0 vibrational frequency 
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FIG. 2: The sketch of the calculated structural and vibrational properties for the top and bridge sites. The reported quantities 
are: dc-o - carbon-oxygen bond length , hco - height of CO molecule above the surface (metal surface - carbon distance), 
dcu-c - Cu-C bond length, Adi2, (resp. Ad23 - average inter- layer spacing between first and second layer), 61,62 - buckling of 
1 st and 2 nd layer, vc-o - stretching mode of C-O molecule where arrows denote distorsions, vs-co - S(urface)- CO vibration. 
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0.0 0.04 
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4.48 4.57 


4.90 


4.65 5.07 
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5.00 



TABLE II: Calculated structural properties of CO adsorbed at the high symmetry sites (t - top, b - bridge) on Cu(211) for 
coverages varying between O = i ML and G = | ML. The reported quantities are: dc-o - carbon-oxygen bond length , hco - 
height of CO molecule above the surface (metal surface - carbon distance), dcu-c - Cu-C bond length, Adi2, Ad23 - change of 
the average inter-layer spacing, 61,62 - buckling of I s * and 2 nd layer, $ - workfunction. The calculated properties are sketched 
in Fig. H 



of 2071 cm- 1 . 

Two molecules in the (3x1) unit cell increase the num- 
ber of possible vibrational modes, adding out-of-phase 
vibrational frequencies. These frequencies are included 
as the second value for the top-top and top-bridge con- 
figurations in Table ITVl We observe two important differ- 
ences between the vibrational frequencies of top-top and 
top-bridge species in the (3x1) surface unit cell: different 
splitting of the vc-o an d vs-CO modes. The in-phase 
symmetric mode vc-o for the top-top configuration is 
calculated to be 2061 cm -1 . The out-of-phase symmet- 
ric mode where both molecules vibrate in opposite di- 
rections is 2007 cm -1 . In the top- bridge configuration 
the in-phase vc-o mode (2057 cm -1 ) is close to the in- 



phase mode of the top-top configuration, but the out- 
of-phase symmetric mode is by almost 100 cm" 1 softer 
(1908 cm- 1 ). 

The experimental EELS study provides values for 
both in-phase and out-of-phase CO stretching modes : 
2088 cm -1 - in-phase and 1881 cm -1 - out-of-phase which 
is already an indication that the adsorption takes place 
in the top-bridge configuration. 

The second difference concerns the splitting of the 
v s-co modes. The in-phase and out-of-phase splitting 
is 2 cm -1 for the top-top configuration and 31 cm -1 for 
the top-bridge configuration. The experimental EELS 
study [33 provides only one S(urface)-CO frequency (341 
cm -1 ) for the high-coverage phase with two symmet- 




■ 

CO on top of Cu 




FIG. 3: STM-images and sphere models of atomic and molecular arrangements during an STM experiment, in which a single 
CO-molecule was pushed along the upper part of a step edge. Notice that a Cu-adatom is located by STM manipulation close 
to the lower step edge (A, a). Upon manipulation along the upper step edge the CO molecule goes to the Cu-adatom (B,b). A 
larger force is necessary to move the CO from the Cu-adatom back to the step edge in accordance with present calculations. 



structural phase 
coverage 
site(s) 


p(3xl)CO 
6=| ML 
top(Cu-adatom) top bridge 


p(2xl)CO 
6=| ML 
top bridge 


p(3xl)2CO 
6=| ML 
top-top top-bridge 


-E a ds,pw9i [eV] 

-Eads.RPBE [eV] 

£ ads ,GGA+u,t/ = L0 eV E e V] 


-1.08 -0.96 -0.98 
-0.84 -0.68 -0.65 
-0.90 -0.72 -0.59 


-0.96 -0.98 
-0.68 -0.65 


-0.94 -0.97 
-0.65 -0.67 



TABLE III: Adsorption energies E a ds of CO on Cu(211) calculated with PW91, RPBE exchange-correlation functionals and 
with the GGA+U method. 



ric modes. In addition to the experimentally observed 
mode at 341 cm -1 we would expect another mode at a 
lower frequency and with appreciably smaller intensity. 
This was not observed in Ref. [33 probably due to its 
small intensity and the large contribution of the specular 
beam in this frequency region. The vs-CO frequency of 
306.7 cm -1 was calculated for CO in off-top site with- 



out phonon bath and the frequency of 343 cm -1 with 
the phonon bath at 100 K 0] • The later frequency is in 
nice agreement with our calculated stretching frequency 
of 346 cm -1 (resp. 342 cm -1 for i ML coverage) and the 
experimentally observed frequency (341 cm -1 ). 
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structural phase 
coverage 
site(s) 


p(3xl)C0 
0=| ML 
top(Cu-adatom) top bridge 


p(2xl)CO 
e = \ ML 
top bridge 


p(3xl)2CO 
0=| ML 
top-top top-bridge 


uc-o [cm -1 ] 
us-co [cm -1 ] 


2071 2049 1925 
346 333 287 


2061 1938 
342 291 


2061,2007 2057,1908 
341,339 334,303 



TABLE IV: Calculated CO frequencies; symmetric stretching mode (yc-o and S(urface)-CO vibration (ys-co)- There are 
two vibrational frequencies, in-phase and out-of-phase, for two molecules in the unit cell. 



D. STM pictures 

Calculated charge density contours have frequently 
been used for comparison with experimental STM images 
and have been able to shed light on complex adsorption 
systems [43|. 

In Fig. 01 we compare typical experimental STM- 
pictures (Fig. 0Ji an d Fig. 0Ji) with calculated charge 
density contours (Fig. 2Jd and Fig. 0^) for both CO- 
configurations discussed for the (3x1) phase, namely top- 
bridge and top-top; the corresponding structure mod- 
els are show in Fig. and f. The calculated top- 
bridge arrangement (Fig. 2b) actually resembles closely 
the experimental STM- picture (Fig. |2Ji). The situation 
is, however, quite different for the case of the top-top- 
arrangement (Fig. QJj and Fig. 0Ji) : Whereas a double 
peak structure is seen in the calculations there is only 
one peak visible in the STM-pictures, whose maximum 
lies in between the two adjacent CO-molecules. 

Furthermore, in most STM pictures the peak appears 
elongated in the direction perpendicular to the steps in 
disagreement with the charge density contour calcula- 
tions. This is alarming, since this top-top-configuration 
is unequivocally established and identified in careful STM 
manipulation and imaging experiments [f|. We can rule 
out experimental reasons for an explanation of this dis- 
crepancy like lateral averaging of the corrugation due to 
final tip width or thermal effects smearing out the double 
peak structure due to vibrations of the molecules. The 
reason for this discrepancy seems more to be caused by 
physical effects which are connected with delicate prob- 
lems in imaging CO-molecules on metal surfaces |4l|: 
It is well known that with metal tips CO is imaged 
as depression (see the isolated CO indicated by an ar- 
row in the STM- image of Fig. E3\). If a CO-molecule 
is deliberately transferred to the tip and CO's at the 
surface are imaged with a CO-terminated tip, then the 
molecules appear as protrusions ■ ^ should be noted 
that charge density calculations for the hypothetical low 
coverage (3xl)-phase yield pronounced maxima at the 
CO-sites. Therefore, isolated CO-molecules imaged with 
a CO-terminated tip resemble more the calculated charge 
density contour. 

The reason for the fact that with clean metal tips 
isolated CO's appear as depressions was recently traced 
back by Nieminen et al. [44j to interference effects be- 
tween different electron paths passing near and through 
the CO, which are canceled when the tip also bears a CO- 



molecule. The STM-picture of the top-top arrangement 
in Fig. 0Jl was taken with a clean metal tip, but also 
with a CO-terminated tip the two adjacent molecules 
appear as a single protrusion. On the other hand, the 
STM-picture of the extended (3xl)-phase (Fig. QJi) was 
highly likely taken with a CO-terminated tip and might 
therefore look more close to the calculated charge density 
contour (Fig. 0^). Nevertheless, in view of the experi- 
mental results for the top-top-configuration (Fig. EJi), 
we cannot strictly exclude an interpretation of the STM- 
picture of Fig. 2Ji, as consisting of top-top CO-pairs sep- 
arated by an unoccupied edge site. On the other hand, in 
view of the present theoretical result, that the top-bridge 
configuration is slightly energetically favored, it appears 
also possible that this configuration is thermodynami- 
cally favored due to the interplay of energy and configu- 
rational entropy effects 0, 0] . More experimental data 
are needed to solve this problem unequivocally. The cal- 
culated differences in work function presented here might 
be of help in this respect. 



V. SUMMARY 

We have used the DFT code VASP to investigate clean 
and CO-covered Cu(211) vicinal surfaces. We consider 
the most probable sites according to experiment: on-top, 
bridge and combinations of these two sites. We show that 
the clean Cu(211) surface tries to smooth its surface and 
relaxes inwards in the first and second layers and then it 
compensates this relaxation by an outward relaxation of 
the third layer. 

CO molecules adsorb perpendicular to the macroscopic 
surface in on-top sites at low coverages up to O = | ML. 
From energy considerations it is non-trivial to distinguish 
which phase is the most favorable for O = | ML. Our 
calculated adsorption energies for the CO molecules sug- 
gest the top-bridge configuration of the CO molecules on 
the Cu(211) at this high coverage. However, the top-top 
configuration is the only one safely established experi- 
mentally. 

From the vibrational properties the top-bridge config- 
uration appears to be thermodynamically favored, but 
more experiments are needed to settle this subtle point. 
Calculations of the CO-adsorption energy on a single 
Cu-adatom adjacent to a step edge show an increase of 
~ 160 meV (RPBE) with respect to CO-adsorption on 
Cu-atoms embedded in the step edges. This is in agree- 
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FIG. 4: Experimental and calculated STM images and corresponding models with two CO atoms in (3x1) cell: (a) experimental 
STM image for top-bridge configuration, (b) calculated STM-image for top-bridge, (c) model for top-bridge configuration, (d) 
experimental STM image for top-top configuration, (e) calculated STM image and (f) model of top-top configuration. The 
colors of the minimal (0 A) and maximal (2.4 A) height of the isosurface are given in (g). The minimal height corresponds to 
the distance of 2.8 A above the surface for the top-bridge and 2.5 A for top-top configuration. 



merit with experimental observations using lateral ma- 
nipulation of the CO. This result is very important for 



vertical manipulation as it shows that CO prefers to ad- 
sorb at low coordination sites as available at the tip apex. 
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